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The electroweak gauge symmetry is allowed to be spontaneously broken by the strongly 
interacting vector-like gauge dynamics. When the gauge coupling of a theory runs slowly 
in a wide range of energy scale, the theory is a candidate for walking technicolor. This 
may open up the possibility that the origin of all masses may be traced back to the gauge 
theory. We use the Schrodinger functional method to see whether the gauge coupling of 
10-flavor QCD "walks" or not. Preliminary result is reported. 
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, 1. Introduction 



The main goal of Large Hadron Collider (LHC) is to confirm the Higgs mechanism 
and to find particle contents and the physics law above the electroweak scale. So far 
many new physics models beyond the standard model have been proposed. Among 
them, Technicolor (TC^ is one of the most attractive candidate^ as it does not 
require elementary scalar particles which cause, so-called, the fine-tuning problem. 
This model is basically a QCD-like, strongly interacting vector-like gauge theory. 
Therefore, lattice gauge theory provides the best way to study this class of model,!^ 
and the predictions can be as precise as those for QCD, in principle. 

The simple, QCD-like TC model, i.e. an SU(3) gauge theory with two or three 
flavors of techniquarks, has been already ruled out by, for instance, the S-parametei^^ 
and the FCNC constraints. However, it has been argued that, if the gauge coupling 
runs very slowly ( "walks" ) in a wide range of energy scale before spontaneous chiral 
symmetry breaking occurs, at least, the FCNC problem may disappear.El Such TC 
models are called walking technicolor (WTC) and several explicit candidates are 
discussed in semi-quantitative manner in Since the dynamics in WTC might 

be completely different from that in QCD and hence the use of the naive scaling 
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in Nc or Nf to estimate various quantities may not work, the S'-parameter must 
be evaluated from the first principlesPI Although really important quantity is the 
anomalous dimension of ^p'tp operator, looking for theories showing the walking 
behavior is a good starting point. Recently many groups started quantitative studies 
using lattice technique to answer the question what gauge theory shows walking 
behavior. In Ref.j^^the running couplings of 8- and 12-flavor QCD are studied on the 
lattice using the Schrodinger functional (SF) scheme.'^ Their conclusion is that while 
8-flavor QCD does not show walking behavior 12-flavor QCD reaches an infrared 
fixed point (IRFP) at ^ 5. In spite of the scheme-dependence of running and its 
value of IRFP, the speculation inferred from Schwinger-Dyson equatiorP^ suggests 
that ^jj^ ~ 5 is not large enough to trigger spontaneous chiral symmetry breaking. 
Although 12-flavor QCD is still an attractive candidate and is open to debate,'^ we 
explore other Nj. In the following, we report the preliminary results on the running 
coupling in 10-flavor QCD. Since the conference, statistics and analysis method are 
changed. The following analysis is based on the increased statistics and a slightly 
different analysis method. 



2. Perturbative analysis 

Before going into the simulation details, let us discuss some results from perturbative 
analysis. In this work, we adopt the P function defined by 



dg\L) 
dL 



= b^g\L) + b2 9\L) + 63 .9" (i) + &4.g'"(i) 



(1) 



where L denotes a length scale. The first two coefficients are scheme-independent, 
and given by 
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Other higher order coefficients are scheme-dependent and are known only in the 
limited schemes and orders. In this section, we analyze the perturbative running in 
the four different schemes/ approximations: i) two- loop (universal), ii) three- loop in 
the MS scheme, iii) four-loop in the MS scheme, iv) three-loop in the Schrodinger 
functional scheme. The perturbative coefficients relevant to the following analysis 
are 
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2-loop universal 








27.74 


9.47 


3.49 


0.52 


3-loop SF 


43.36 


23.75 


15.52 


9.45 


5.18 


2.43 


0.47 


3-loop MS 




159.92 


18.40 


9.60 


5.46 


2.70 


0.50 


4-loop MS 






19.47 


10.24 


5.91 


2.81 


0.50 



Tab. 1 The IRFP from perturbative analysis. 



where the coefRcients and Cg depend on the spatial boundary condition of the 
SF used in calculations, i.e 9. Those for 6 = tt/^ and C2 for 6 = arc known as 

^9=^/5 ^ 125563 + 0.039863 x Nf, (6) 

£3=''/^ = {4=""^^^ + 1.197(10) + 0.140(6) X Nf- 0.0330(2) x 7V|, (7) 

c^=° = 1.25563 + 0.022504 x Nf, (8) 

but C3=° is not. Therefore, the case iv) is studied with 6=71/5. Notice that in our 
numerical simulation 0=0 and thus, rigorously speaking, the example iv) is not 
applied to our numerical result. 

The perturbative infrared fixed point (IRFP) are numerically solved and sum- 
marized in Tab.[T] As seen from Tab.[TJ with Nf > 8 the fixed point value is, to some 
extent, stable against the change of schemes/approximations. It is interesting that 
the IRFP of 12-flavor QCD in the SF scheme is consistent with the non-perturbative 
calculation by Ref.l^Now looking at the perturbative IRFP at Nf = 10, it appears 
to be stable at g"^ ^ 10. Furthermore, according to the analysis based on Schwinger- 
Dyson equation, there is an argument that chiral symmetry breaking occurs at 
around ~ 4^V(3 C;2(i?)) = ^^Pln summary, the perturbative analysis suggests 
that 10-flavor QCD is the most attractive candidate for WTC. 

3. Simulation parameters and setup 

We employ the Schrodinger functional methocP to calculate the running coupling 
constant. Unimproved Wilson fermion action and the standard plaquette gauge 
action without any boundary counter terms are used. The parameter of the spatial 
boundary condition for fermions, 6, is set to zero. The bare gauge coupling /3 = G/g^ 
is explored in the range of 4.4-24.0. In this analysis, we report the results obtained 
from [L/a)^ = 6"^, S"' and 12^ lattices. The calculation on 16^ lattice is in progress. 
The numerical simulation is carried out on several architectures including GPGPU 
and PC cluster. The standard HMC algorithm is used with some improvements 
in the solver part like the mixed precision algorithm. So far, we have accumulated 
5,000 to 200,000 trajectories depending on (/?, L/a). 

Since the Wilson type fermion explicitly violates chiral symmetry, the critical 
value of K, has to be tuned to the massless limit. We performed this tuning for 
every pair of {(i, L/a). At around /3 4.4, for massless fermions we encounter a 
(probably first order) phase transition independently of L/a, where the plaquette 
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value suddenly jumps to a smaller value. Since this bulk phase transition is inferred 
to be lattice artifact, whenever this happens we discard the configurations. Thus 
the position of the critical /3 (~ 4.4) sets the lower limit on /3 at which simulations 
make sense. 



4. Preliminary results 

Figure [T] shows the Schrodingcr functional coupling calculated on the lattices, where 
ffo/ffsF plotted as a function of the bare coupling. The solid curve is the result 
(and statistical error) of the fit to 

gg ^ 1 - flL/aa go 

gSF^(50' ^/«) 1 + Pl-L/a X gl + E^=2 Oi/a^n X " ' 

where Pi,L/a is the L/a-dependent coefficient and is found, by perturbative calcu- 
lation, to be 

r 0.4477f 07831 for L/a = 6 
Pi,L/a =< 0.4624813408 for L/a = 8 . (10) 
i 0.4756888260 for L/a = 12 

We optimize the degree of polynomial in the denominator of ([9|) by monitoring 
xVdof, and take TV = 5 for L/a = 6 and TV = 4 for L/a = 8, 12. 

Since we do not implement any 0{a) improvements, large scaling violation is 
expected to exist. One promising prescription to improve discretization errors has 
been proposed in Rcf.'^ Let us parameterize the lattice artifact in the step scaling 

by 

^^{u,s,L/a) ^^^^^^ 

!i Ti \ 1 + 5^'^\s,L/a)u ' ^(2), Ti \ 2 . 

5{u, s, L/a) = ^ ' = 5^^' {s, L/a)u'^ H , (11) 
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Fig. 2. 6 as a function of u. S^'^'l (u, s, L/a) is determined as a coefficient of the quadratic term in 
u. (solid and dashed curves are obtained from different fit ranges.) 



where u = gspi^), (t{u,s) = glpisL) and I]o(w, s, L/a) is g^^j {g^,sL/a) at g^ 
satisfying g^'l^ {gl,L/a) = u. 

The coefficient S^^''{s,L/a) in eq. (fTTj) is given by 

^(1) [s, L/a) = (p^,sL/a - bi ln(sL/a)) - (p,^L/a - bi ln(L/a)) . (12) 

Dividing the lattice data T,o{u,s,L/a) by the factor (1 + S^^^^ {s , L / a)) improves 
the 0{u) discretization error and hence 6(u,s,L/a) starts from O(u^) as already 
indicated in eq. (fTTj) . 

(t(u, s) has perturbative expansion, 

(t(u, s) = U + SqU^ + SlU^ + S2U'^ + S3U^ + • • ■ , (13) 

so = 61 ln(s), si = ln(s) {h^ ln(s) + 62) , (14) 

52 = ln(s) (^61^ ln2(s) + ^6162 ln(s) + b^j , (15) 

53 = ln(s) l^i^ ln3(s) + ^bi%2 In'(s) + ln(s) (^86163 + ^62') + , (16) 

where fe^'s are the coefficients of the /3- function introduced in secO Since we know 
the first two coefficients bi and 62 and hence cr(u, s) can be numerically determined 
to 0(u"^), with such a the 0{u^) term in S{u, s, L/a) is attributed to discretization 
error. The coefficient of term, (^'■^'(s, L/a), is then obtained by fitting S{u, s, L/a) 
to a quadratic function of u. This fit must be done in the small- coupling region where 
the perturbative series is reliable. Since at this moment we have only a limited 
number of data points in such a region, the fit range is forced to extend to u 2.5. 
Figure [2] shows the u dependence of S{u, s, L/a) and the fit results. The obtained 
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Fig. 3. Discrete beta functions for [s, L/a)=(4/3, 6), (2, 6) and (3/2, 8) from left to right. Two 
colored bands are the results with statistical error, obtained by the two-loop improvement with 
eq. I I17I1 (red) and eq. I I18I1 (blue), respectively. The black solid (dashed) curve is the correspond- 
ing discrete beta function obtained by integrating perturbative two-loop universal (three-loop SF 
scheme) beta function. 



coefficients are 

r 0.0062(9) for s = 4/3, L/a = 6 
S'^^'> (s, L/a) = < 0.0138(26) for s = 2, i/a = 6 , (17) 
[ 0.0070(24) for s = 3/2, L/a = 8 

r 0.0053(8) for s = 4/3, L/a = 6 
(5(2)(.s, L/a) = <^ 0.0122(21) for s = 2, i/a = 6 . (18) 
i 0.0062(19) for s = 3/2, L/a = 8 

In the fit, two different fit ranges, u g [0, 2.02] and [0, 2.50], are applied to examine 
the fit range dependence of the result. Eqs. ([T7|) and eq. ([T8|) correspond to the 
former and the latter fit range, respectively. Using 5^'^\s,L/a) thus extracted, we 
define the improved lattice data by 

= i4-Mi)r f/^I'^l'^ilt ri ^ 2 - (19) 

1 -t- o^^iys, L/a) u + b^'^>[s, L/a) 
To see the running in detail, we introduce the discrete beta functioUjli^ 

i?'^* {u,s, L/a) = \ - - . (20) 

Figure [3] shows the 1/u dependence of B^^^{u,s,L/a). As usual /3 function, large 
negative value of i?'*^* means rapid increase with length scale of the coupling, and 
^lat fjippjjig the sign indicate the existence of IRFP. As seen from the figure, in 
either pairs of (s, L/a) the discrete beta function approaches to zero from below 
when 1/u decreases from 0.8 to 0.2, and this happens independently of the choice 
of S'-^K This means that at around u=:1.25 the running starts to slow down, i.e. 
"walk" . In order to confirm that this behavior persists even in the continuum limit, 
simulations on a larger lattice is necessary. 



5. Summary 

The running coupling constant of 10-fiavor QCD is studied. The perturbative anal- 
ysis suggests that this theory is extremely interesting. The preliminary result ob- 
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tained without continuum limit seems to indicate the walking behavior. In order 
to draw definite conclusions, we clearly need larger lattices to take the continuum 
limit. Such calculations are in progress. 

A part of numerical simulations is performed on Hitachi SRllOOO and the 
IBM System Blue Gene Solution at High Energy Accelerator Research Organiza- 
tion (KEK) under a support of its Large Scale Simulation Program (No. 09-05), 
on GCOE (Quest for Fundamental Principles in the Universe) cluster system at 
Nagoya University and on the INSAM (Institute for Numerical Simulations and Ap- 
plied Mathematics) GPU cluster at Hiroshima University. This work is supported in 
part by the Grant-in- Aid for Scientific Research of the Japanese Ministry of Educa- 
tion, Cuhure, Sports, Science and Technology (Nos. 20105001, 20105002, 20105005, 
21684013, 20540261 and 20740139), and by US DOE grant #DE-FG02-92ER40699. 
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